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Abstract. We show that non-uniqueness of the Leray-Hopf solutions of the Navier- 
Stokes equation on the hyperbolic plane H 2 observed in |CCj is a consequence of 
the Hodge decomposition. We show that this phenomenon docs not occur on H n 
whenever n > 3. We also describe the corresponding general Hamiltonian setting 
of hydrodynamics on complete Ricmannian manifolds, which includes the hyperbolic 
setting. 

Introduction 

Consider the initial value problem for the Navier-Stokes equations on a complete 
n-dimensional Riemannian manifold M 

(1) d t v + Vyv — Lv = — gradp, divv = 

(2) v(0,x) = v o (x). 

The symbol V denotes the covariant derivative and L = A— 2r where A is the Laplacian 
on vector fields and r is the Ricci curvature of M. Dropping the linear term Lv from 
the first equation in §1§ leads to the Euler equations of hydrodynamics 

(3) d t v + Vyv = — gradp, divv = 0. 

Most of the work on well-posedness of the Navier-Stokes equations has focused on the 
cases where M is either a domain in IR n or the flat n-torus T n . In fundamental contribu- 
tions J. Leray and E. Hopf established existence of an important class of weak solutions 
described as those divergence-free vector fields v in L°°([0, oo), L 2 ) n L 2 ([0, oo), if 1 ) 
which solve the Navier-Stokes equations in the sense of distributions and satisfy 

(4) \\ v (t)\\ 2 L 2 + A f ||Def v(s)\\ 2 L 2ds < \\v \\ 2 L 2 and lim \\v(t) - v \\l^ = 

Jo 

for any < t < oo and where Def v = |(Vf + Vf T ) is the so-called deformation tensor. 
When n = 2 using interpolation inequalities and energy estimates it is possible to show 
that the Leray-Hopf solutions are unique and regular but the problem is in general 
open for n — 3, see e.g. |CF] or |MB] . 
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There have also been studies on curved spaces, which with few exceptions have 
been confined to compact manifolds (possibly with boundary), see e.g. |Taj and the 
references therein. In a recent paper Chan and Czubak [CCj studied the Navier- 
Stokes equation on the hyperbolic plane H 2 and more general non-compact manifolds 
of negative curvature. In particular, using the results of Anderson |Anj and Sullivan 
|Su] on the Dirichlet problem at inifnity, they showed that in the former case the 
Cauchy problem (JTJ- (J2]) admits non- unique Leray-Hopf solutions. 

Our goal in this note is to provide a direct formulation of the non-uniqueness of the 
Leray-Hopf solutions on M 2 which turns out to rely on the specific form of the Hodge 
decomposition for 1-forms (or vector fields) in this case. We also show that no such 
phenomenon can occur in the hyperbolic space H n with n > 3. As a by-product, we 
describe the corresponding Hamiltonian setting of the Euler equations on complete 
Riemannian manifolds (in particular, hyperbolic spaces). 

We point out that this type of non-uniqueness cannot be found in the Euler equa- 
tions. Furthermore, it is of a different nature than the examples constructed e.g., by 
Shnirelman |Sh] or De Lellis and Szekelyhidi |DSj . On the other hand, it is similar 
to non-uniqueness of solutions of the Navier-Stokes equations defined in unbounded 
domains of the higher- dimensional Euclidean space, cf. Heywood |Hej . 

1. Stationary harmonic solutions of the Euler equations 
Our main result is summarized in the following theorem. 
Theorem 1.1. 

(i) There exists an infinite- dimensional space of stationary L 2 harmonic solutions 
of the Euler equations on M 2 . 

(ii) There are no stationary L 2 harmonic solutions of the Euler equations on M n 
for any n > 2. 

Proof. Recall the Hamiltonian formulation of the Euler equations j3} on a complete 
Riemannian manifold M, see e.g. |AKj . Consider the Lie algebra Q Teg = Vect /J (M) of 
(sufficiently smooth) divergence-free vector fields on M with finite L 2 norm. Its dual 
space Q* cg has a natural description as the quotient space VL l L2 / d£l° L 2 of the L 2 1-forms 
modulo (the L 2 closure of) the exact 1-forms on M. Namely, the pairing between cosets 
\fi] G Vt l L2 /dVt Q L 2 of 1-forms /3 G Vt l L2 and vector fields w G Vect M (M) is given by 

(\J3],w) := / (i w (3) dfi, 

J M 

where l w is the contraction of a differential form with a vector field w, and /i is the 
Riemannian volume form on M. 

Let A : g rcg — > Q* cg denote the inertia operator defined by the Riemannian metric. 
The operator A assigns to a vector field v G Vect M (M) the coset [v ] of the corresponding 
1-form •& via the pairing given by the metric. The coset is defined as the 1-form up 
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to adding differentials of the L 2 functions on M. Thus, in the Hamiltonian framework 
the Euler equation reads 



where [v^] G Q} L2 /dVt° L2 and L v is the Lie derivative in the direction of the vector field 
v. 

The space of the L? 1-forms on a complete manifold M admits the Hodge- 
Kodaira decomposition 

n\ 2 =~^®8n 2 ~;®H 1 L 2, 

where the first two summands denote the L 2 closures of the images of the operators d 
and S, while is the space of the L 2 harmonic 1-forms on M. Therefore, we have a 
natural representation of the dual space 



d; eg = sQ 2 L2 ®nl 2 . 

It turns out that the summand of the harmonic forms in the above representation 
corresponds to steady solutions of the Euler equation. Namely, one has the following 
proposition. 

Proposition 1.2. Each harmonic 1-form on a complete manifold M which belongs to 
L 2 fl L A defines a steady solution of the Euler equation Q on M. 

Proof of Proposition Let a be a bounded L 2 harmonic 1-form on M. Let v a denote 
the divergence-free vector field corresponding to a, i.e., v° a = a. Since the 1-form a is 
harmonic, using Cartan's formula gives 

^-a = -L Va a = -L Va da - di Va a = -di Va a. 

We claim that t Va a G f2° 2 and consequently dt Va a G dQ° L2 - Indeed, by the definition 
of the vector field v a we have 

\\ L v a at\\ 2 L 2 = (a(v a )) 2 dfi= \\a\\li, 
Jm 

which is finite by assumption. It follows that the 1-form dt Va a must correspond to 
the zero coset in the quotient space $j* eg = VL l L2 /dVt Q L2) which in turn implies that 
■^a = G 0* eg . The latter means that the 1-form a defines a steady solution of the 
Euler equation, which proves the proposition. □ 

If M is compact then the space of harmonic 1-forms is always finite-dimensional 
(and isomorphic to the deRham cohomology group H 1 (M)). According to a well- 
known result of Dodziuk [Doj . the hyperbolic space HP carries no L 2 harmonic fc-forms 
except for k = n/2, in which case it is infinite-dimensional. Therefore, there can be no 
L 2 harmonic stationary solutions of the Euler equations on HP for any n > 2, which 
proves part (ii) of the theorem. 

To prove part (i) we note that for n = 2 the space of harmonic 1-forms on H 2 is 
infinite-dimensional. Moreover, it allows for the following construction. Consider the 
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subspace S C of 1-forms which are differentials of bounded harmonic functions 
whose differentials are in L 2 

S = {(/<£> | $ is harmonic on H 2 and #6i 2 }. 

It turns out that the subspace S is already infinite-dimensional. Indeed, let us consider 
the Poincare model of HI 2 , i.e., the unit disk D with the hyperbolic metric ( , )h, 
which we denote by B>h- It is conformally equivalent to the standard unit disk with 
the Euclidean metric ( , ) e , denoted by D e . Bounded harmonic functions on can 
be obtained by solving the Dirichlet problem on D e , i.e., by constructing harmonic 
functions $ on D with boundary values <p prescribed on <9D. First, the 1-form ci$ is 
clearly harmonic: 

Arf$ = d6d$ = rfA$ = 0. 

Secondly, observe that 



5;, ) 



(d$,d<S>) h dfi h 
(ci$,ci$) e ci/i e 



det(g ij ){d<f>, d$) e det(^) dfi e 



\d*\\ht 



and 



If/* I 



(d$,d$)ldfi h = / det (g v )(d^ 1 d^} e det(g ij ) d\i e 



|2\2, 



(d$,d<5>) 2 e dfi e (z) < / rf$) 2 dfi e = 



|4 

li 4 (D e )' 

where det(gij) = 1/(1 — \z\ 2 ) 2 is the determinant of the hyperbolic metric. 

Furthermore, for sufficiently smooth boundary values </? G C 1+a (dB>) there is a uni- 
form upper bound for its harmonic extension inside the disk: 

\d${x)\ < C\\(p\\ci+*{m) 

for any igD and < a < 1, and some positive constant C, see e.g. |GT] . This implies 
that (for sufficiently smooth tp) the 1-forms d& define an infinite-dimensional subspace 
S of harmonic forms in L 2 flL 4 , which satisfy assumptions of the proposition above. It 
follows that they define an infinite-dimensional space of stationary solutions of the Euler 
equations on the hyperbolic plane H 2 . This completes the proof of Theorem 11.11 □ 



2. Non-unique Leray-Hopf solutions of the Navier-Stokes equations 

Using the fact that suitably rescaled steady solutions of the Euler equations also 
solve the Navier-Stokes system the authors in |CC] obtained a type of ill-posedness 
result for the Leray-Hopf solutions in the hyperbolic setting. 

Theorem 2.1 QCCj ). Given a vector field v e = (ci$)" on H 2 there exist infinitely many 
real-valued functions f(t) for which v ns = f(t)v e is a weak solution of the Navier-Stokes 
equations with decreasing energy (i.e., satisfying the Leray-Hopf conditions). 
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An immediate consequence of this result and Theorem 11.11 is the following 

Corollary 2.2. There exist infinitely many weak Leray-Hopf solutions to the Navier- 
Stokes equation on M 2 . There are no non-unique Leray-Hopf harmonic solutions to the 
Navier-Stokes equation on HP with n > 3 arising from the above construction. 

Remark 2.3. The phenomenon of nonuniqueness of solutions to the Navier-Stokes 
equation in unbounded domains Q C K n , n > 3, of higher- dimensional Euclidean 
spaces is of similar nature, see |Hej . Indeed, that construction is based on the existence 
of a harmonic function with gradient in L 2 and appropriate boundary conditions in 
such domains. The Green function $(x) = G(a, x) centered at a point a outside 
of Q has the decay like G(a, x) ~ |x| 2 ~ n as x — > oo, so that |<i$(x)| ~ Ixl 1 "™ and 
hence IgM^x)! 2 ~ |x| 2 ~ 2n . Thus, for n > 3 the 1-forms c?$ belong to L 2 ft L A on fl. 
The corresponding divergence- free vector fields (d§Y provide examples of stationary 
Eulerian solutions in Q (with nontrivial boundary conditions) and can be used to 
construct time-dependent weak solutions v ns = f(t)(d$Y to the Navier-Stokes equation 
in Q, as in Theorem 12. II 



3. Appendix 

To make this note self-contained we provide here some details of the construction 
of the weak solutions given in |CCj . It will be convenient to rewrite the Navier-Stokes 
equations ([T|) in the language of differential forms 

(5) d t v b + V v v b - Av b + 2r(w b ) = -dp, 5v° = 

where Sv 9 = — div v and Av b = d5v^ + 5dv ^ is the Laplace-deRham operator on 1-forms. 

Let v be the vector field v ns = on H 2 as in Theorem 12. II Since the 1-form 

d$> is harmonic one only needs to compute the covariant derivative term and the Ricci 
term: 

V^l = \f 2 {t) dm 2 and 2r(v b ns ) = -2f{t)d*. 

Direct computation, taking into account the fact that for H 2 we have r = — 1, shows 
that both terms can be absorbed by the pressure term, so that the pair (i^p), where 
p := (2/(t) - /'(*))$ - l/2/ 2 (t)|rf$| 2 satisfies the equations ©. 

Finally, a quick inspection shows that any differentiable function f(t) satisfying 

f 2 {t)+A f f 2 (s)ds<f 2 (0) 
Jo 

yields a vector field v ns which satisfies the remaining conditions in (jlj) required of a 
Leray-Hopf solution. 
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